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This review is dedicated to the memory of two 
pioneers, Teijo Aberg and Dick Deslattes, who led, 
each in his own way, generations of researchers to a 
deeper and more comprehensive understanding of the 
Physics ofX Rays. 

1. Introduction 

A wavelength conveniently close to the repeat 
distance of many crystals, reasonably low absorption in 
low- and medium-Z elements, and favourable physical 



Present address: X-FAB Semiconductor Foundries AG, Erfrirt, 
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and chemical properties facilitating the construction of 
x-ray tubes conspired to make the Cu Ka x rays the 
most widely used radiation in diffraction and scattering 
experiments. These wide-employed applications, and 
the intriguing asymmetry of the Cu Ka and Kj8 line- 
shapes [1,2], along with those of all l>d transition 
elements [3-8], led, in turn, to a century of extensive 
spectrometric studies of the Cu Ka and Kj8 spectra. In 
spite of these extensive studied, recent studies reveal 
that surprises still lurk under the skewed Kai2 and 
overlapping Kj8i 3 lines, and the related multi-electron- 
ic satellite (S) and hypersatellite (HS) spectra. 

The asymmetric lineshape of the copper emission 
lines were attributed in the past to a number of differ- 
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ent processes: Kondo-like interaction of the conduc- 
tion electrons with the core-holes [1,2,9], final-state 
interactions between the core holes and the incomplete 
3d shell [6], 2pl3d shell electrostatic exchange inter- 
action [7] and, most importantly, shake-up and shake- 
off of electrons from the 3/ shells [8,10]. The last 
process, in particular, received in the past strong exper- 
imental support [10-12]. Our measurements of both 
Kai2 and K^Si 3 diagram spectra to be described below 
[13-15], carried out in conjunction with relativistic 
multiconfigurational Dirac-Fock (RMCDF) ab initio 
calculations, also supports this interpretation, and 
provide accurate numerical estimate of the magnitude 
of the contribution from these effects to the lineshape. 
Moreover, we have carried out recently excitation- 
energy-dependent measurements of the spectra near the 
excitation energy threshold for these shake-off effects. 
Our results lend further support to the assignment of the 
asymmetry of the Cu diagram lines to 3d spectator-hole 
transitions (in the following, underlining denotes hole 
states). 

To support a detailed and accurate shape analysis, 
high resolution measurements of the spectra are 
required. This, in turn, requires a detailed study of 
the various factors affecting the resolution of the 
measuring system, and their optimization for a given 
measurement. Thus, in preparation for the Cu Kai 2 and 
Kj8i3 measurements, we have carried out a careful 
study and characterization of our measuring system 
[16]. This allowed to minimize the spectral distortions 
due to finite resolution. Also, our use of a calibrated sil- 
icon analyzer allowed an uncertainty on the order of 
10"^ (Kj8) or better (Ka) in the absolute energy scale. 
The use of the same well characterized experimental 
setup under well defined resolution conditions for both 
the Ka and Kj8 spectra provided higher confidence in 
the conclusions drawn from the data. We have also car- 
ried out ab initio relativistic Dirac-Fock calculations to 
help determine the structure underlying the lines. These 
assumed that the spectral lineshapes can be accounted 
for by satellites resulting from 3/ spectator hole transi- 
tions in addition to the nominal single electron diagram 
transitions . We have found a specific set of transitions 
which accounts extremely well for the Ka, and to a 
somewhat lesser extent, also for the Kj8 lineshapes 
[13-15]. 

The theoretical understanding of the diagram lines, 
even when considering only the dominant single- 
electron-transition components, becomes more compli- 
cated when the near-threshold region is considered. In 



the "isothermal" region, high above threshold, the 
excitation and de-excitation processes of inner shell 
electrons have been extensively studied [17], and can 
be well described theoretically by the prevailing 
Sudden- or Frozen-Core approximation, where the 
ejected electron is removed immediately, the atom's 
shell structure is kept frozen in its ground state config- 
uration, and the electrons treated as independent, non- 
interacting and uncorrected entities [18]. The excita- 
tion and de-excitation processes are treated in this 
regime as two independent and consecutive processes. 
By contrast, in the "adiabatic" regime [19], near the 
excitation threshold, effects like inter- and intra-shell 
electronic correlations, atomic shell relaxation and the 
changing interaction between the slow-moving ejected 
electron and the relaxing atom, become progressively 
more important, even when single-cXcctxon transitions 
are involved [18,20,21]. Near threshold, the excitation 
process, and the de-excitation by photon (xrays) or 
electron (Auger) ejection can no longer be considered 
to be independent processes; they become increasingly 
simultaneous, mutually interacting, and merge into a 
single complex process. In this regime the Sudden and 
Frozen-core approximations are no longer valid in 
principle. Thus, the study of the near threshold region 
is very important for elucidating electronic correlations 
and the interplay between excitation, relaxation and 
emission in atoms. 

To study these effects, we have used tunable 
synchrotron radiation to excite the diagram lines at 
near-threshold energies and follow the evolution of the 
spectral shape and intensity from below threshold to 
saturation [22,23]. These measurements, reported 
below for the first time, indicate a threshold for the 
asymmetry of the lines at -10 eV above the K edge. 
Below this threshold the Cu Ka and Kj8 lines are found 
to have symmetric Lorentzian shapes. The asymmetry 
was found to grow with excitation energy and to 
saturate -100 eV above threshold. These results strong- 
ly support the assignment of the asymmetry to 3d spec- 
tator transitions, the threshold of which is expected to 
occur for Cu at ^8 eV to 10 eV above the K edge, 
according to the Z+1 approximation [24]. 

Spectator-hole transitions, like the 3d spectator-hole 
transitions contributing to the diagram lines, involving 
more than a single electron or vacancy, go beyond the 
simple single electron, frozen atom description of the 
excitation and emission processes. They require taking 
into account correlations between electrons in the 
same, or in different shells [18]. Such multielectronic 
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effects can best be studied in satellite (2/?-spectator 
transitions, denoted Ka3 4, originating in ls2p ^'2p^ 
transitions) and hypersatellite (l^-spectator transitions, 
denoted K^ori,25 originating in 1^^ ^ ls2p transitions) 
spectra, which , on the one hand, originate in transitions 
involving two vacancies in their initial and final states, 
and, on the other hand, are removed fi'om, and thus not 
masked by the strong diagram lines, as are the 3d spec- 
tator lines straddling the Ka and Kj8 diagram transition 
lines. The study of satellite and hypersatellite spectra, 
particularly near threshold, was hampered until very 
recently by the lack of suitable excitation sources 
which are tunable, narrow-band, and intense enough to 
allow studying these weak transitions. Pioneered by 
Deslattes and coworkers [25], and fiirther developed by 
others, synchrotron-based beamlines equipped with 
efficient, high resolution fluorescence spectrometers 
have become available over the last few years [26,27]. 
Such facilities allow detailed studies of this kind, as 
well as inelastic scattering experiments, to be carried 
out [22,25,28]. We will demonstrate here the richness 
and sophistication achievable in studies using this 
powerfiil technique by discussing briefly two of our 
recent studies: the CuKa34 satellites [22,23], and the 
CuK^ai^2 hypersatellites [29,30]. 

We have studied the variations in the Cu Ka34 
emission spectrum (ranging fi-om 8060 eV to 8100 eV) 
as a fimction of the photoexcitation energy from 
threshold at ===10000 eV up to the saturation of the 
satellites' intensity at 1 1 200 eV. The results confirm 
the pure shake-off nature of the spectrum, as predicted 
by theory. They also reveal two distinct regimes in the 
spectral evolution. In the first, up to ===50 eV above 
threshold, the spectral shape as well as the intensity 
undergo a rapid and complicated variation with excita- 
tion energy. In the second regime, fi'om -50 eV to 
==1000 eV above threshold, no variation of the spectral 
shape is observed, and only the overall intensity 
increases monotonically to saturation at the range's 
upper end. The saturated spectral shape, though not the 
fast-varying one near threshold, is found to be in good 
agreement with ab initio RMCDF calculations, which 
allows to identify the various spectral features with 
specific transitions. The Thomas Model [31,32], devel- 
oped specifically for the cross-section evolution in the 
"adiabatic" regime near threshold employing time- 
dependent perturbation theory, though successfiil in the 
case of valence electrons in low-Z atoms [33], does not 
agree well with our experimental results, indicating a 
strong need for more detailed theoretical work in this 
almost unexplored regime. 



Hypersatellites (HS), are of particular interest for a 
number of reasons. First, since the initial state involves 
two electrons in the same shell, the spectra should 
allow studying intra-sheW correlations. Also, since 
even for medium-Z atoms the K shell is already rather 
strongly relativistic the HS allow studying the effects of 
relativity on electronic correlations in atoms. 
Moreover, the K^ai HS line originates in the ^Sq ^^Pi 
spin flip transition, which is dipole-forbidden in the 
pure LS coupling scheme and fully allowed only in the 
jj coupling scheme [24,34]. Thus, the intensity ratio 
R=I( K^ai)//(K^a2) is strongly dependent on the 
degree of intermediacy of the coupling [35], rendering 
it the most sensitive (and almost only) method for 
quantitatively studying the coupling variation across 
the periodic table fi*om the LS coupling limit at low Z to 
the jj coupling limit at high Z. K^a^ 2 HS spectra are also 
unique in allowing to study the Breit interaction, the 
most elusive and least studied of all atomic inter- 
actions. Because of a near cancellation of the Coulomb 
terms, the contribution of the Breit interaction to the 
splitting of the K^ 0:1,2 HS lines , and to their shift fi'om 
the diagram lines, can reach as high as 20 % at high Z, 
rather than the small 1 % contribution in the case of 
diagram lines [36,17]. Finally, HS are the "diagram" 
transitions of, and thus allow studying hollow atoms 
[37], i.e., atoms where a whole inner shell is empty, 
while the outer shells are occupied [38]. The formation 
mechanism of such atoms and their properties are 
of fundamental interest to basic atomic physics 
[37,40]. They are also of great importance for studies of 
atoms very far fi-om equilibrium and of ultra-fast 
dynamics in atoms, with possible wide-ranging appli- 
cations in physics, chemistry, biology and materials 
science [41], as well as a possible way towards the 
realization of hard-x-ray lasers [41,42]. 

Using the same technique employed for the Cu 
satellites, we have measured the first clean, resolved 
Cu K^ai2 hypersatellite spectrum [29,30]. The line 
splitting and shift from the diagram line agree well with 
our, and previous, RMCDF calculations. Nevertheless, 
these calculation slightly overestimate the measured 
intensity ratio R. A typical shake-off behaviour with 
intensity increasing smoothly from threshold is found, 
in agreement with shake theory predictions and previ- 
ous satellite measurements. The Thomas Model does 
not agree well with the measured intensity evolution, as 
also found for satellites. The evolution of the spec- 
trum's intensity from threshold shows an unexpectedly 
large saturation range, ==60 % of the threshold energy. 
Taken together with the results, mentioned above. 
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obtained for the M and 2p spectator transitions, a dis- 
tinct trend emerges in the length of the saturation range 
with orbital quantum number /. To the best of our 
knowledge, such a trend has not hitherto been 
addressed by theory. 

In the following we discuss the experimental and 
calculational methods, the results obtained for the dia- 
gram lines, and, briefly also the satellite and hyper- 
satellite results. 

2, Experiment and Calculations 

2.1 Single and Double Crystal Laboratory-Based 
Spectrometry 

The laboratory-based diagram line measurements 
were carried out using a conventional sealed x-ray tube 
source, operated at 40 kV and 30 mA. For accurate 
wavelength measurements a single-crystal spectrome- 
ter was employed, while for the lineshape measure- 
ments a double-crystal spectrometer was used. Both 
spectrometers were optimized using the detailed theo- 
retical and experimental study of Hartwig et al. [16]. 

The single crystal spectrometer allowed absolute 
angle measurements with a mean total angle dividing 
error of 0.12" arc-sec and an angular step size of 0.06" 
[43]. Using the WAS09 silicon crystal [44], previously 
calibrated at the German institute of standards, PTB, 
this spectrometer allowed, in principle, wavelength 
measurements to an accuracy of AA/A = 3 x 1 0"^ 
[44,45]. For the CuKa spectrum the symmetrical (444) 
reflection was used and for the CuK^S spectrum- the 
asymmetrical (553) reflection. The incident beam was 
defined by a 390 mm long collimator, with horizontal 
slits 0.42 mm high and 0.03 mm (CuKa) and 0.10 mm 
(CuKjS) wide. Under these conditions the measured 
intensity distribution is well approximated by a convo- 
lution of the true spectral distribution and the (approx- 
imated) instrumental function. The true distribution 
was obtained in this case by deconvolution of the 
measured spectrum. 

The double crystal spectrometer is a commercial 
DTS spectrometer [46]. Since the zero angle of 
this device could not be determined accurately, 
the energy scale was fixed by defining the angular 
position of the maximum of the CuKai peak to be A^ = 
0.154 059 292 nm [44]. For the double crystal spectro- 
meter it was possible to make the instrumental function 
very narrow, so that no correction for finite resolution 
was necessary, and only a correction for absorption was 
apphed to the data [16]. Here symmetric Si(333) 



reflections were used in both crystals for both spectra. 
A 560 mm long collimator was employed with a 
horizontal slit width of 1.0 mm (focus side) and 10 mm 
(sample side) for both spectra, and vertical slit heights 
of 0.4 mm (CuKa) and 0.8 mm (CuK^. 

Further details of the experimental setup of the dia- 
gram line measurements are given in Refs. [13,15]. 

2.2 Photoexcited Synchrotron-Based Spectrometry 

The measurements in the threshold region of the 
diagram lines, and all the satellite and hypersatellite 
spectra, were carried out at beamlines X2 1 and X25 at 
NSLS, Brookhaven National Laboratory. The experi- 
mental setup is shown schematically in Fig. 1. In the 
X25 beamline, for example, the source is a hybrid 
27-pole wiggler. The radiation fi-om the source is 
focussed by a toroidal mirror, followed by a 2-bounce 
Si (220) or a Si (1 1 1) monochromator. This provides an 
incident energy resolution of 2 eV to 6 eV, and a flux of 
-lO^^photons/s in a spot size of < 0.7 mm^ at the 
sample. The incident intensity is monitored by an ion- 
ization chamber. The fluorescence radiation fi^om the 
sample is measured by a Johann type spectrometer, 
having a Rowland circle of 1 m diameter on a hori- 
zontal plane, and a spherically bent 3 in diameter Si 
analyzer. For the diagram Kai2 (at -8050 eV), the 
satellite Ka;34 (at ^8080 eV), and the hypersatellite 
K^ai^2 (at -8330 eV) spectra the Si(444) reflection was 
used. For the diagram K^ 3 (at -8900 eV) the Si(553) 
reflection was used. For the energies studied the 
analyzers have high Bragg angles, ranging from 70° to 
80"^, which, in conjunction with the small focus of the 
exciting radiation, provide for a resolution of 2.8 eV or 
lower for these spectra. A liquid nitrogen cooled intrin- 
sic Germanium detector was used to measure the ana- 
lyzer-focussed radiation. Further details, especially 
on the beamline and spectrometer, measurement 
procedures and data treatment are given in [22,23,26]. 

2.3 Ab Initio Transition Calculations 

The calculations employed in the data interpretation 
in all studies were done ab initio, using the relativistic 
multiconfigurational Dirac-Fock (RMCDF) package 
GRASP [47], with supplementary code written in- 
house. Following previous work [21,48] the initial and 
final state wavefiinctions were generated in separate, 
independent runs. The energies of the individual transi- 
tions are then obtain by subtracting the appropriate 
level energies, as calculated in the initial- and in the 



78 



Volume 1 09, Number 1 , January-February 2004 

Journal of Research of the National Institute of Standards and Technology 




Two — bounce 
mono. 



Spherical 
anolyzer 



Sannple N;^ 



Detector 



Fig.l. Schematic view of tlie synclirotron-based setup for measuring 
near-thresliold pliotoexcited emission spectra. The incident exciting 
beam is monochromatized by a double-bounced silicon crystal 
monochromator. The emitted radiation is recorder by a Johann-type 
spectrometer using a spherically-bent analyzer in near-back-reflec- 
tion geometry. 



fmal-state runs. This procedure was found to yield 
accurate energies for the CuKa34 satellites even for 
excitation energies within 50 eV of the threshold [22]. 
The calculation of relative transition probabilities 
within each multiplet requires the wavefunctions of the 
initial and final states to be orthogonal. Since the initial 
and final states are generated here in separate runs, this 
condition is not fulfilled. Thus, as in previous studies 
[13,22,49], configuration interaction calculations were 
carried out to obtain the various transition probabilities 
using once the initial state orbitals and again those of 
the final state. All the significant transition probabilities 
in the two sets agreed with each other to within ± 10 %. 
We used, therefore, the line strengths calculated from 
the initial state wavefiinctions. 

2.4 Phenomenological and Theoretical Spectral 
Fits 

To obtain an accurate analytical, albeit phenomeno- 
logical, representation of the measured data the dia- 



gram lines were phenomeno logically fitted by a sum of 
Lorentzians, using non-linear least-squares methods. 
The number of Lorentzians in the sum was increased 
until a good % ^ value, which did not decrease further 
upon increasing the number of Lorentzians, was 
obtained. Such a representation is particularly impor- 
tant for the diagram Ka and Kj3 spectra, which are used 
in diffraction and scattering experiments. 

For interpreting the structure of the spectra, the 
"stick diagram" multiplets obtained in the calculations 
above were also fitted to the measured data. Here, each 
transition line was replaced by a single Lorentzian, the 
integrated area of which was equal to the relative inten- 
sity of the corresponding "stick" within the multiplet. 
Usually, all Lorentzians representing lines within a 
given multiplet were assigned a common width, 
although in some cases (particularly for the diagram 
spectra interpretation) two independent widths were 
used: one for transitions associated with the a^ or j3i 
lines, and the other for transitions associated with the 
a^ or ^83 lines. A spectrum was typically fitted by a sum 
of several multiplets, where for each multiplet a single 
intensity parameter was refined in the fit, representing 
its relative intensity contribution to the given spectrum. 
Further details are given in [13,15,25]. 



3. Results and Discussion 

3.1 The Kai,2 and KjSj^ "Diagram" Spectra 

3.1.1 Analytical Representation of the Spectra 

The Ka and Kj3 , measured using the double crystal 
spectrometer as discussed above, are shown in Figs. 2 
and 3. The basic parameters derived from the measured 
spectra [15] are summarized in Table 1. A fiill analysis 
of this data is given in Ref [15]. Here we mention only 
the ±3x10"^ (Kai,2) and ±4x10"^ (KA,3) relative 
accuracy of the absolute wavelengths, the highest 
to date for these lines. The removal of the finite 
resolution effects by optimizing the spectrometer to 
a level of a negligible resolution width (in the case 
the double crystal spectrometer) or by the decon- 
volution of the data by a well-determined resolution 
fiinction (in the case of the single-crystal spectrometer) 
yields accurate intrinsic lineshapes, which allow 
carrying out a detailed theoretical analysis, as discussed 
below. 
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Fig. 2. Double-crystal spectrometer measured and smoothed CuKa3 4 spectrum (points) and fitted Lorentzians. 
The fit parameters are given in Table 2. The lower panel shows the fit residuals. The thin lines denote the ±2<j 
values of the data, where a is the standard deviation. 
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Fig. 3. Same as Fig. 2 but for CuKft 3. 



Table 1, The peak wavelengths on an absolute scale, A , energies, E, full widths at half maximum, W, and asym- 
metry indices, a^, of the measured diagram lines. Measurement uncertainties are given in brackets. Wavelength 
to energy conversion is done using 806 554.477(32) (eVm)~^ [75]. Note the 10"^ level of uncertainty of the 
present wavelength results. 



Line 


A(nm) 


E(cV) 


ff(eV) 


Ow 


CuKcfi 


0.154 059 29 (5) 


8047.83 (1) 


2.29 (2) 


1.07 


CuKa2 


0.154 442 74 (5) 


8027.85 (1) 


3.34(6) 


1.36 


CuKjSi3 


0.139 223 4 (6) 


8905.42 (4) 


5.92 
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For analytic representation purposes, the Ka 
spectrum was fitted by four Lorentzian, as shown in 
Fig. 2. The Lorentzians are defined by their position E^p 
width Wy and relative intensity /^, where / = 1 ,2 denotes 
Kai2, respectively andy is the index within each line. 
The fit parameter values are given in Table 2. The fit 
residuals (lower panel in Fig. 2) are small and non- 
systematically scattered around zero. Moreover, they 
are within the ±2(7 limits of the data , where c denotes 
the standard deviation of the data points due to count- 
ing statistics. This residual distribution indicates 
an excellent fit, also demonstrated by the weighted 
7?-factor of 0.7 % obtained in the fit. 

The double-crystal-spectrometer-measured Kj3 
spectrum is shown in Fig. 3. As noted earlier [50] four 
underlying lines are discernible. These were denoted 
as the two diagram lines j3i and ^83 and two satellites; 
j3 ' on the low energy side, and j3^^on the high energy 
side. A fifth line was proposed by Bremer and S0rum 
[12], but no determination of its position, intensity and 
width was carried out. Hayasi [51] also observed this 
feature and determined its energy without a detailed fit. 
Bremer and S0rum also observe a very sharp dip just 
off the peak, on the low energy side, where only a 
smooth shoulder was observed in earlier measure- 
ments. This was assigned to the spin-orbit splitting of 



the two diagram lines. No such dip was observed in any 
of the other published line shapes, nor in our measure- 
ments. Finally, although the position of the Kj3 ' feature 
as measured by Salem et al. [52] is in good agreement 
with theory, its intensity disagrees with both the 
exchange interaction [52,53] and plasmon oscillation 
[54] theories. 

To address these issues, we have carried out a series 
of fits with four and five Lorentzians [13]. We found 
that indeed the use of five Lorentzians yields a signifi- 
cant improvement in the fit residuals and R factors over 
the four-Lorentzian one [13]. The parameters defining 
the five-Lorentzian representation of the lineshape, as 
obtained from the fit, are listed in Table 2. Although the 
weighted 7?-factor of this fit, 1 %, is somewhat higher 
than that of the Ka spectrum, the fit is still very good, 
and , as can be seen in the lower panel of Fig. 3, the 
residuals are within the ±2(7 limits of the measured 
points. The association of the Lorentzians j3a to ^^ with 
the phenomenologically determined features j3i, ^83,^3 \ 
and j3^^of the spectrum is discussed in Ref [13]. 

Further details of the various quantities derived 
directly from the raw Ka and Kj3 spectra, the phenom- 
enological resolution of these spectra into sums of 
Lorentzians, and their comparison with previous 
measurements are given in Refs. [13,15,16]. 



Table 2. Positions Ey, widths Wy, amplitudes ly and integrated intensities l^^^ obtained for tlie individual Lorentzians from a fit of each meas- 
ured diagram line by the sum of the minimal number of Lorentzians which provides a reasonable fit to the data. Weighted 7?-factors, R^ , obtained 
for the fits are also shown. 



Line 



Comp. 



^/,(eV) 



Wy{cV) 



^integ 



CuKai2 



«21 
«22 



8047.837(2) 


2.285(3) 


0.957(2) 


8045.367(22) 


3.358(27) 


0.090(1) 


8027.993(5) 


2.666(7) 


0.334(1) 


8026.504(14) 


3.571(23) 


0.111(1) 



0.579 
0.080 
0.236 
0.105 



^w(%) 



0.7 



CuKj8i3 



8905.532(2) 


3.52(1) 


0.757(3) 


8903.109(10) 


3.52(1) 


0.388(2) 


8908.462(20) 


3.55(3) 


0.171(2) 


8897.387(50) 


8.08(8) 


0.068(1) 


8911.393(57) 


5.31(8) 


0.055(2) 



0.485 
0.248 
0.110 
0.100 
0.055 



^w(%) 



1.0 
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3.1.2 Comparison With ^6 Initio Theory 

To elucidate the structure underlaying the diagram 
lines , we carried out RMCDF calculations for transi- 
tions originating in several different initial states. 
In additional to the diagram transitions Is ->2p we 
considered also the two-hole transitions ls3l ^2p3l 
where l = s,p,d, and the most probable three-hole 
transition \s3d^ ^2p3d^, The resulting "stick" dia- 
grams of these multiplets are given in Figs. 4 and 5 for 
the Ka and Kj3 spectra, respectively along with the 
measured spectra and the individual Lorentzians' 
positions in the phenomenological fits described above. 
One can observe in these plots that all multiplets 
considered are well-aligned with the spectra (within the 



< 2 eV computational accuracy [13]), and thus can, in 
principle contribute to the measured spectrum. Since 
the lifetime widths of the transitions are larger than the 
interline distances, each of these multiplets will result 
in a smooth lineshape of overlapping contributions and 
it is not possible to identify by eye the individual con- 
tributions of the various multiplets. Thus, only general 
conclusions can be drawn from eye inspection. 
For example, the strongest lines of the \s3p -^ 2p3p 
multiplet are in a region where the measured intensity 
is small, and hence the contribution of this multiplet to 
the line shape can not be large. A visual inspection of 
Fig. 5 also reveals that the 3d spectator hole multiplet 
has strong lines in the vicinity of the P \ P\, and j3^2 
features, as pointed out by LaVilla [11]. 
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To obtain more quantitative conclusions, we have 
carried out a series of computer fits, where each "stick" 
in a given calculated multiplet was represented by a 
Lorentzian, whose integrated intensity (area under the 
curve) is equal to the height of the corresponding 
"stick". The widths of all Lorentzians within a multiplet 
were constrained to be equal. A series of fits with dif- 
ferent combinations of the diagram lines with one or 
more multihole multiplets were carried out, refining the 
width, overall intensity and position of each multiplet. 
The details of the fits are given in Tables III and VII in 
Ref [13] and will not be reproduced here. The R factors 
and Goodness-of-Fit indicators presented there lead 
clearly to the conclusion that the inclusion of the 3d 
spectator transition, ls3d -> (2p, 3p)3d, in the fit along 
with the diagram ones. Is -^ (2p, 3p), is mandatory for 
obtaining a good fit for both of the measured Ka and 



Kj3 spectra. No indication of contributions from the 3^ 
spectator transitions was found. This is not surprising in 
view of the fast depopulation of the 3^ vacancies by 
Coster-Kronig transitions. Similarly no indication of 
contributions from the 3p spectator transitions were 
found for Kj3, though a very small contribution, -0.5 % 
of the total intensity to the Ka spectrum may be possi- 
ble. The best fit, that including the diagram and the 3d 
spectator transitions is shown in Figs. 6 and 7 for the 
two spectra, along with the residuals. Note that all 
residuals of the Ka are within ±2(7, and are randomly 
distributed. For the Kj3, some deviations of the residu- 
als from this limits occur, due to the P ' (-8897 eV) and 
the j3''(-8908 eV) satellites [11]. However, the small 
residuals and the good weighted R factors, 3 % for Ka 
and 5.3 % for Kj3, clearly support the conclusion 
that the only appreciable contribution to the lines 
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fit residuals are shown in the lower panel. 



besides the diagram transitions is from the 3d spectator 
hole transition, in agreement with earher resuhs [10,7] 
based on non-relativistic HF calculations. 

Finally, we note that the best fit, shown in Figs. 6 and 
7, yields a rather high contribution of the 3d spectator 
transition to the total intensity of the spectrum: 26 % 



to 30 % in both spectra. This demonstrates clearly 
that multielectronic effects play a more important 
role than hitherto assumed in determining the spectral 
shape of even the strong diagram lines. Similar or even 
stronger influence can be expected in the weaker 
atomic spectra. 
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3.1.3 The Near-Edge Evolution of the Spectra 

The interpretation above for the structure underlying 
the Ka and Kj3 spectra entails different thresholds for 
the symmetric (diagram) and the asymmetric (3d 
spectator) contributions to the lineshapes. The sym- 
metric part is due to a single-electron transition, origi- 
nating in a single Is vacancy Thus, its threshold is at 
the binding energy of a single Is electron in a neutral 
Cu atom. This is the conventional K edge, at 8979 eV 
[55]. The asymmetric part, due to the 3d spectator 
transition \s3d ^ (2p, 3p)3d originates in a double 



ls3d vacancy. Using the Z+1 approximation [24] and 
the ^8 eV binding energy of the 3d electrons in Zn [55], 
the threshold for the 3d spectator transition, the KM 
edge, should be about 8 eV to 10 eV above the K edge. 
An immediate implication of this conclusion is that 
below this KM edge, but above the K edge, the specta- 
tor transition is not excited and the Cu Ka and Kj3 line- 
shapes should be symmetric and well-represented by a 
single Lorentzian for each transition. Two points need 
be noted. Since the 3d electron binding energy rises fast 
with Z from Cu up, the actual KM threshold may occur 
at an energy higher than 8 eV above the K edge. 
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Moreover, as shake-off may dominate the creation of 
the initial \s3d vacancy state, and since the cross- 
section for a shake-off process rises gradually from 
zero at threshold, the contribution of the 3d spectator 
transition also rises gradually from threshold. 
Immediately above threshold the contributions from 
the 3d spectator transitions are still small. Thus, the 
asymmetric features of the lineshape may not become 
discemibly until a somewhat higher excitation energy 
is reached. 



We have measured the evolution of CuKa and Kj8 
spectra as the photoexciting energy varies from the K 
edge and up, using beamline X21 at NSLS, 
Brookhaven National Laboratory, and the experimental 
setup described above. The results for the Ka^ line are 
shown in Fig. 8, and those for the overlapping Kj8i3 
complex in Fig. 9. The additional structure, which caus- 
es the lineshapes' deviation from simple Lorentzians 
and eventually makes them asymmetric, can be clearly 
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Fig. 8. The evolution of the CuKa^ Une with photoexciting energy near threshold from an almost symmetric 
single-Lorentzian lineshape (soUd line) to an increasingly asymmetric lineshape. The asymmetry results from 
the openning and subsequent growth of the 3d spectator hole transition channel starting from 10 eV to 15 eV 
above the K edge [E (K edge) = 8979 eV]. 



87 



Volume 109, Number 1, January-February 2004 

Journal of Research of the National Institute of Standards and Technology 



300 



250 



'""'^"'=9000 eV 
=9050 eV 
X-ray tube, 40 kV 
Single-electron 
(calculated) 



8 200 

CO 



CO 

o 



on 



50 



00 



50 






8885 8890 8895 8900 8905 8910 8915 8920 

ENERGY (eV) 

Fig. 9. Same as Fig. 8, but for CuKjSi 3. Here the lineshape below and at threshold of the 'id spectator hole 
transitions is that of two overlapping Lorentzians (solid line) as expected of the two overlapping diagram 
transitions KjSj and Kj33. 



seen to increase from an excitation energy of Eg^citation^ 
8890 eV and up. As observed in Fig. 10, particularly in 
the very small fit residuals in the lower panel of the 
figure, the Kj3i 3 spectrum excited at ^'excitation = 8890 eV 
can indeed be excellently fitted by two Lorentzians 
only, representing the two one-electron Kj3i and Kj33 
diagram transitions: \s -^3p 



Careful examina- 
tion of the growth of the 3d spectator features in both 



1/2, 3/2- 



spectra with increasing excitation energy reveals that 
the multielectronic features grow monotonically, and 
reach their maximal intensity about 100 eV to 150 eV 
above the KM threshold. This range is in good 
agreement with recent measurements of a 250 eV to 
300 eV growth range of the GeKj3'^^ satellites, which 
originates in similar 31 {l = s,p,d) spectator transitions 
[56]. 
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Further details of our measurements of the near-edge 
evolution of the spectral features of the Ka and Kj8 
spectra, with special emphasis on the evolution of 
features associated with spectator vacancy transitions, 
were published elsewhere [57]. 

3.1.4 Monolithic Silicon Spectrometer Measurements 

We finish this brief review of our measurements of 
the CuKa and Kj8 spectra by a short discussion of a 
novel, though by now not so new, monolithic crystal 



device allowing absolute wavelength determinations at 
1 0"^ level at every laboratory without the use of a pre- 
viously calibrated crystal. 

The Monolithic Double Crystal Spectrometer 
(MDCS) [58-60] is shown in Fig. 11. It is an x-ray 
optical element cut from a single block of a perfect 
crystal. Two sets of Bragg planes of the monolith play 
the roles of the two separate crystals of a conventional 
double-crystal spectrometer (DCS) [61]. The require- 
ment of satisfying two Bragg conditions simultaneous- 
ly determines the transmitted wavelength, given the 
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Fig. 11. The monolithic double crystal spectrometer (MDCS). a is 
the angle of rotation and {hy ky ly) and (/?2 ^2 ^ ^^"^ the two partici- 
pating Bragg planes. The beam trajectory used for aligning the 
rotation axis normal to the (hy ky ly) plane, and the trajectory of the 
measured beam are also shown. 



angle between the two reflecting planes. Energy scan- 
ning is done by rotating the MDCS by an angle a about 
the normal to the first set of Bragg planes. This leaves 
the angle of incidence (Bragg angle) on the first set of 
planes unchanged, but varies the angle of incidence on 
the second set, and consequently also the wavelength 
transmitted by the MDCS. This method of wavelength 
scanning was pioneered in conventional double crystal 
spectrometry by Deslattes [62] and employed since in 
numerous other studies [63]. Nevertheless, the first 
detailed theoretical study of the MDCS was published 
only recently [60]. 

For an x ray to be transmitted through the MDCS at 
a given angle of rotation a, its wavelength has to be 
[60] 

Xia)/d, =2sin^,cosa/[(7(/^'+^'+0/(/?>^f+/f) 
-cos^o)'+sin'^,cos'a]''' (1) 

here dy is the spacing of the first of the two sets of 
planes, (hy, ky, ly) and (^2, ^2? h)^ and Pq is the dihedral 
angle between these two planes. Two important points 
should be noted. First, the wavelength transmitted is 
directly related in Eq. (1) to the lattice constant 

^0 =^i^J^ +]<^ + '^ . Since particular commercial 



brands of single crystal silicon boules are now 
produced routinely to part-per-million tolerances for 
the lattice constant, cutting an MDCS from such a 
boule allows absolute wavelength determination of an 
x-ray line to 10"^ with moderate effort and simple 
equipment available in standard x-ray laboratories [58] 
without having to resort to calibrated crystal slabs 
obtained from (few and often busy) standards labora- 
tories. For example, the particular MDCS used 
to measure the spectrum in Fig. 12 below yields 
A (CuKai)/ao= 0.283 665 49(24), i.e., an uncertainty of 
0.85 X 10^ [58]. Second, in a conventional DCS the 
transmitted wavelength A ©^ sin a and, consequently, 
the wavelength change upon rotation of the second 
crystal is linear in the angle of rotation: 
AA oc Aa. For the MDCS, however, the transmitted 
wavelength according to Eq. (1) is approximately 
A oc cos a so that for small rotation angles the change 
in the wavelength is quadratic in the angle of rotation: 
AA oc (Aa)^. A significantly larger rotation is required, 
therefore, for the MDCS than for the DCS to obtain the 
same wavelength change. This provides considerably 
higher angular resolution and dispersion for the MDCS. 
The CuKai line, measured using a silicon MDCS with 
(111) and (206) planes, is shown in Fig. 12. The 
quadratic dependence of the energy on the rotation 
angle is clearly observed. This quadratic dependence of 
AA on a for the MDCS has several other implications. 
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rotations a of the MDCS in Fig. 11, cut from a perfect silicon 
crystal using the (HI) and the (206) planes. 
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First, the wavelength transmitted depends on the 
absolute value of a only, not on its sign. Thus, the same 
spectrum is measured twice; once for the negative and 
once for the positive senses of rotation. This symmetry 
relative to a = allows to determine accurately the 
zero point of the rotation from the measured data, for 
example, as the mid-point between the same peak 
measured for +a and -a rotations. This double meas- 
urement also allows the identification of features 
appearing in one spectrum only as spurious, so that it 
can be excluded from the consequent analysis. These 
advantages come at a cost of a limited scan range, and 
the existence of a high limit for the measurable 
wavelength, X^ , corresponding to the zero angle, a = 0, 
position of the MDCS: 



XJd.^l sin % l[{4{^+1^+ll)l{h\+k\+ll) 

-cos I3,y+ sin' I3,f\ (2) 

This limit has to be kept in mind in choosing the 
Bragg planes to be used in an MDCS designed for 
scanning a given energy range. To obtain the highest 
dispersion and resolution, the planes have to be chosen 
such that the longest wavelength of the range is as close 
as possible to, but still smaller than, Xq. Further details 
can be found in Refs. [58,59], which also discuss in 
detail the sources and magnitudes of the experimental 
errors of the MDCS. A detailed theoretical analysis of 
the MDCS is given in Ref [60]. 

3.2 The K(X3,4 Satellite Spectrum 

Using the same technique employed for the near- 
threshold region of the diagram Ka and Kj8 spectra, we 
have measured the spectral evolution of the CuKa34 
satellite complex from threshold up. The measurements 
were carried out at beamline X25, NSLS, Brookhaven 
National laboratory. The full results of that study have 
been published recently [23,23], and will be reviewed 
here only briefly. 

A two regime behaviour was found in the evolution: 
in the first ^50 eV above threshold, a fast variation of 
both the shape and intensity of the spectrum is found. 
Above that, only the intensity keeps growing monoton- 
ically with excitation energy, without any changes 
occurring in the spectral shape, until the intensity 
saturates at ^1 keV above the threshold. We now 
discuss each regime separately. 

In Fig. 13 we plot the background-subtracted 
spectrum measured in the high-energy regime. We also 
plot the RMCDF ab initio calculated transition lines 
corresponding to the 2p spectator transitions 



ls2p ^'2p', which were suggested already in 1927 
[64,65] to be the originators of the ^cc^^4 spectrum. The 
overall alignment of the calculated lines with the meas- 
ured spectrum and, in particular, its various features is 
good. The four main features, a\ a^, a^, and a 3 can 
clearly be identified with the ^P^ ^^P^, ^P2^^P2, 
^Pi ^ ^D2, and ^Pi^^P2 transitions. We have also 
measured the energy threshold for exciting each of the 
four main features and found them to be 10017 eV, 
10011 eV, 10027 eV and 10015 eV, for the a\ a^, a^, 
and a'^ lines, respectively [23]. The ab initio "stick" 
diagram multiplet of the 2p spectator transitions was 
fitted to the measured spectrum, using a Lorentzian to 
represent each transition (i.e., each "stick"). Two fits 
are shown in Fig. 13. Fit A employed a single width 
common to all the multiplet lines, and only this width, 
a single intensity scale factor, and a small computation- 
al shift (~1 eV to 2 eV) between the energy scales of 
the measurement and the calculation were refined in the 
fit. Fit B allowed the widths and relative intensities of 
the individual "stick" to vary independently, and only 
the individual line positions within the multiplet were 
held fixed. The agreement of even the more restricted 
fit, Fit A, with the measured spectrum is good, indicat- 
ing that no other multiplets contribute to the lineshape. 
Indeed, attempts to include contributions from the 
calculated 2s spectator transitions 1^2^ -^2s2p invari- 
ably reduced their intensity to zero [22]. This is in line 
with the ~5 fold lower shake probability calculated for 
a 2s electron, as compared to a 2p one, to accompany a 
\s vacancy production \66'\. The very strong Coster- 
Kronig transition 1^2^ -^ Is2p31, which depopulates 
the 2s spectator state very fast, fiirther reduces any 
possible 2s spectator contributions to the spectrum 
[13,67]. Finally, we wish to point out that the fact that 
the as-calculated ab initio multiplet (Fit A) requires 
adjustments, albeit small, of the individual lines' 
integrated intensities (via their widths) to achieve a 
near-perfect fit (Fit B) indicates that while our calcula- 
tion captures the essentials of the spectrum, further 
effects, such as slightly less-than-fiill relaxation and/or 
final- and initial-state correlations, may need to be 
included in the ab initio calculations to achieve a good 
agreement with the measurements. 

Within 50 eV of the threshold, the shape of the 
spectrum varies rapidly. Several of the spectra meas- 
ured in this regime are shown in Fig. 14. Each spectrum 
is normalized to its maximal intensity, so that the 
weakest spectrum, at ^g^^itation = 1 00 1 e V, is about 
30-fold less intense than the strongest one, at ^excitation = 
10250 eV. This intensity decrease, and the high back- 
ground due to the Ka^ tail, effectively prevented 
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Fig. 13. The measured CuKa34 satellite spectrum (points), with its fit by the ab-initio RMCDF calculations. The fits are 
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the second panel are marked by the LS notation of their major constituents. The calculated 1^2^ -^Islp transitions are 
shown in the third panel. The fits reveal that these transitions do not contribute to the lineshape to any detectable level. 



meaningful satellite spectrum extraction at lower 
^excitation values. Note that for these reasons, and the 
finite energy spread in the exciting beam, remnants of 
some of the emission lines are still observed at the 
lowest ^excitation spectrum, measured at the nominal 
threshold of 10 010 eV. We have used a 4-Lorentzian 
phenomenological fit to allow following the 



variation of each underlying line, since the weaker lines 
could not be refined with any confidence level. For 
obvious reasons we have not attempted to reproduce 
the complicated variation of the lines' intensities with 
incident energy, observed in this figure, by ah initio 
calculation. This would require including inter-shell 
electronic correlations and other near-threshold effects 
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Fig. 14. The measured CuKa34 spectra in the near threshold adiabatic region, for several excitation energies, 
each normalized to its maximum by the scale factor listed. Fits to a sum of four Lorentzians, and the individual 
Lorentzians, are shown in solid lines. Note the roughly constant widths and positions of the lines, and the 
increase of the a^ line from zero below threshold to domination at ~10 080 eV. 



in the calculations. We hope that the availability of our 
measurements in this "adiabatic" regime will spur 
theoretical developments which will eventually allow 
ab initio calculations to be carried out even in this 
highly complex regime. 



Further details, in particular of the variation of the 
satellites' intensity with exciting energy from threshold 
to saturation, are discussed in Refs. [22,23]. Similar 
measurements for satellites of Ge were published 
recently by Stememann et al. [56] 
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3.3 The CuK^'oi^ Hypersatellite Spectrum 

The CuK^ai2 hypersatellite (HS) measurements 
were carried out by beamline X25 with the same setup 
as those of the CuKa34 satellites discussed above. A 
full discussion is given in Refs. [29,30], and only a 
short summary is presented here. 

The HS spectrum measured at ^'excitation = 20 keV is 
shown in Fig. 15. The most outstanding feature of the 



spectrum is the inverted Kai/Ko^ intensity ratio, -0.29, 
as compared to the well-known value of ~2 observed 
for the corresponding diagram lines. This results from 
the fact that in the LS coupling, prevailing for low-Z 
atoms, the K^a^ transition is spin-flip forbidden. It is 
frilly allowed only in the jj coupling, prevailing for 
high-Z atoms. This makes this ratio the most sensitive 
measure for the intermediacy of the coupling at a given 
Z [68-70]. The value we measured for Cu, 0.29(2), is 
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Fig. 15. The measured CuK ai2 hypersatellite spectrum (points) with its fit (solid line) by the RMCDF 
calculated "stick diagram" in the third panel. The fit residuals are shown hi the second panel. Note the good fit 
when QED corrections are included in the calculations. The spectrum calculated without these corrections 
(fourth panel) is upshifted relative to the measured spectrum. It can not be fitted to the measured spectrum even 
when downshifted, due to its too large line splitting (dashed line). This demonstrates the large contribution of 
QED effects to the hypersatellite spectra. 
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slightly overestimated by the best RMCDF calculations 
[68], including ours [30], albeit by no more than two 
standard deviations. 

We have carried out ab initio RMCDF calculations 
of the HS spectrum, both including and excluding QED 
corrections. The closed 3<i shell of Cu and the empty K 
shell of the initial state results in a sparse calculated 
"stick" diagram, shown in the third and fourth panels of 
Fig. 15. Representing each line in this diagram by a 
Voigt function with a fixed Gaussian resolution width 
of 2.8 eV, we fitted the calculated spectra to the meas- 
ured one. The resultant fit for the calculation which 
includes QED is shown in a solid line in Fig. 15. The 
residuals are shown in the second panel of the figure. 
They are almost all within the ±2 cr limits of the meas- 
urements. Thus, within the measurements' statistics, the 
intrinsic lineshape is well described by a single 
Lorentzian, indicating no contamination by higher- 
order multivacancy transitions, such as the 3d spectator 
transitions, which were shown above to contribute 
26 % to 30 % to the total intensity of the corresponding 
Ka^ 2 diagram lines. This is not surprising in view of the 
very low probability for creating the two-K-hole initial 
state of the HS transition, 10^ relative to that of creat- 
ing a single K hole, and the fact that the spectator states 
require one additional hole, i.e., a three-vacancy initial 
state. The vanishingly low probability of exciting such 
states is responsible, then, for the high purity of the 
photoexcited HS spectra. The fit to the calculation 
which neglects the QED corrections is shown in dashed 
line in the first panel of Fig. 15. As can be seen, the 
incorrect splitting does not allow to obtain a good fit to 
the measured spectrum even when the calculated spec- 
trum is shifted to lower energies by the large amount of 
^1 eV required. This demonstrates the exceptionally 
large contribution of QED effects to the HS, and hence 
the importance of the HS spectra for studying these 
effects [68-70]. Comparing the measured spectrum 
with more detailed calculations, which include the var- 
ious components of the QED corrections one at a time 
should elucidate the relative importance of each of 
these contributions. Such calculations are however out- 
side the scope of our (mostly experimental) study. 

The evolution of the peak intensity of the K^a^j line 
with the incident energy, ^'excitatiom is shown in the upper 
panel of Fig. 16. While it seems to saturate at 
^ 23 keV, it should be noted that this curve is propor- 
tional to the total probability of obtaining the initial 
two-K-hole state in a neutral atom, not to the probabil- 
ity of obtaining the second K hole per directly ionized 
first K hole created. The last quantity, which is the one 



most often cited in the literature, has also been evaluat- 
ed [30], and is shown in the lower panel. A fine scan of 
the threshold region is shown in the inset of the 
upper panel, yielding £'^,3hoid= (18.352 ± 0.015) keV, 
in excellent agreement with our Z+1 approximated [24] 
^fcrtlhoid = 18.345 keV and, to a lesser extent, with our 
RMCDF-calculatedCeshoid= 18.378 keV. As can be 
observed, even on this magnified scale the intensity 
rises from threshold smoothly, as expected of a shake- 
off effect, and no jumps, characteristic of shake-w/? 
transitions, are discernible. This is in line with shake 
theory's prediction of an increasing shake-off/shake-up 
ratio with increasing Z and decreasing principal quan- 
tum number, n, of the shaken electron's shell [67,71]. 
Trends in recent DF calculations for noble gases 
[71-73] interpolated to the present case of Z= 29 and 
n= \ predict a contribution of less than 1 % from 
shake-up to the total shake probability. The error bars 
of the data in Fig. 16 yield an experimental upper limit 
of <3 % on shake-up contributions to the HS spectrum 
at threshold in our case. 

The intensity of the CuKa34 satellites, discussed 
above [22], which originate in an n = 2 spectator hole 
transition, were found to saturate at ~1 keV above 
threshold, i.e., -10% of the threshold energy. Here, 
where the spectrum originates in an n = 1 spectator hole 
transition, saturation is not reached even at the highest 
energy measured, 25 keV. Extrapolation [30] indicates 
that saturation is reached only around -30 keV, making 
the saturation range here a surprisingly large -11 keV 
or almost 60 % of the threshold energy. These percent- 
age values, when considered with the 3 % saturation 
range of the /2 = 3 spectator transitions in Ge [56], and 
the 1 % to 2 % range inferred above (from our near- 
threshold measurements of the lineshapes of the Cu 
Kai 2 and Kj3i 3 diagram lines) for the 3<i spectator hole 
transitions [57] indicate a very fast increase in the 
saturation range with decreasing n of the spectator hole. 
We also note that the currently prevailing Thomas 
model [31,32,74] developed to account for the intensi- 
ty variation in the adiabatic regime, does not describe 
well the intensity variation of our measured data in 
Fig. 16: when the threshold is fixed at the measured 
value the shape of the fitted curve varies greatly from 
the measured one, and when the shape is optimized, let- 
ting the threshold vary, the refined threshold is -1 keV 
lower than the measured one. A better understanding of 
the trends in the saturation range's dependence on n and 
/ of the spectator hole, and the variation of the intensi- 
ty with incident energy, must await better theoretical 
calculations in the adiabatic regime. 
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Fig. 16. The intensity variation of tlie CuKa 2 line with incident energy (upper panel). The inset shows a finer 
measurement of the threshold region, yielding the indicated threshold energy. The lower panel shows the 
relative intensity variation of the hypersatellite spectrum (points), along fits of the Thomas model discussed in 
the text (lines). Note that saturation is not reached even at 25 keV, i.e., -6.7 keV above threshold. 



4, Conclusion 

The measurements presented above demonstrate the 
importance of multielectronic transitions in determin- 
ing the structure of the emission hnes of Cu, including 
the Kai2 and Kj3i3 diagram ones. In conjunction with 
detailed ab initio relativistic Dirac-Fock calculations 



we have shown that such transitions can contribute as 
much as -30 % of the intensity of the diagram lines and 
are responsible for their well-known asymmetric 
lineshapes. Inner-shell spectator transitions produce 
satellite and hypersatellite lines which have been 
explored as a function of excitation energy using photo- 
excitation by monochromatized synchrotron radiation. 
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These studies reveal intriguing new effects such as line- 
shape variations near threshold, unexpectedly long 
(and, to date, unexplained) saturation ranges strongly 
dependent upon n and / of the spectator hole, sensitive 
dependence of particular features on the coupling 
scheme and intermediacy etc. These effects, and, in 
particular, the near-threshold behaviour, are still far 
from being understood theoretically It seems therefore 
that in spite of an extensive, century-long activity in 
this field, much still remains to be elucidated. We 
look forward to new discoveries and the emergence of 
deeper understandings in this field in the fiiture. 
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